Molecular spintronics: Coherent spin transfer in coupled quantum dots 
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Time-resolved Faraday rotation has recently demonstrated coherent transfer of electron spin be- 
tween quantum dots coupled by conjugated molecules. Using a transfer Hamiltonian ansatz for the 
coupled quantum dots, we calculate the Faraday rotation signal as a function of the probe frequency 
in a pump-probe setup using neutral quantum dots. Additionally, we study the signal of one spin- 
polarized excess electron in the coupled dots. We show that, in both cases, the Faraday rotation 
angle is determined by the spin transfer probabilities and the Heisenberg spin exchange energy. By 
comparison of our results with experimental data, we find that the transfer matrix element for elec- 
trons in the conduction band is of order 0.08 eV and the spinFaraday rotation signal as a function 
of the probe frequency in a pump-probe setup transfer probabilities are of order 10%. 
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PACS numbers: 78.67.Hc,85.65.-|-h,85.75.-d 



I. INTRODUCTION 

The past years have evidenced rapid experimental 
progress in the field of spintronicsi^ Coherent trans- 
port of electron spins in semiconductors has been demon- 
strated over several micrometers ^2. nourishing hopes that 
the electron spin may be used as carrier of information 
similar to the electron charge. Such applications of the 
spin degree of freedom for classical or quantum informa- 
tion processing^ require control of the electron spin not 
only in extended systems such as two-dimensional elec- 
tron gases (2DEG's), but rather also for spins localized 
in quantum dots (QD's). 

Recently, coherent transfer of electron spin has been 
observed between QD's with different radii ~ 1.7 nm 
(QD A) and tb — 3.5 nm (QD B) coupled by a benzene 
ring.^ The different QD radii give rise to different quan- 
tum size levels for electrons and holes in the two species of 
QD's. The resulting difference in exciton energies allows 
one to pump and probe selectively the spin polarization 
for QD's of species A and B. The main result of Ref. 
is that an electron spin polarization created by optical 
pumping in QD B is transferred 'instantaneously' to QD 
A. The efficiency of this transfer mechanism is of order 
10% at low temperatures T < 50 K and increases to ap- 
proximately 20% for T > 100 K. The observed shift of 
the exciton energies to lower values compared to isolated 
QD's is also consistent with a coherent delocalization of 
the electron or hole over the system formed by the QD's 
and the bridging molecule. 

The purpose of this paper is to show that a two-site 
Hamiltonian with a transfer term captures some of the 
essential experimental features. We aim at calculating 
the dependence of the experimentally observed Faraday 
rotation (FR) signal as a function of probe energy on 
microscopic parameters such as spin transfer probabili- 
ties. The FR angle is proportional to the difference in 
refractive indices for circularly polarized light which 
is determined by the difference of the dielectric response 



functions. We calculate the dielectric response functions 
of coupled QD's and derive an analytical expression for 
the FR angle in terms of electron transfer probabilities 
and Heisenberg exchange splittings. The experimental 
data provide strong evidence that the spin transfer is 
mediated by the 7r-conjugated moleculei^ We do not aim 
to describe this transfer mechanism microscopically, but 
consider the transfer matrix elements for electrons and 
holes as parameters of the Hamiltonian. 

For CdSe QD's with radii rA and rs, the single-particle 
level spacing for electrons and holes is large compared 
to the temperatures T < 200 K explored experimentally. 
This allows us to restrict our attention to the lowest or- 
bital levels in the conduction and valence band of both 
QD's. A possible admixing of higher orbital levels caused 
by the Coulomb interaction is determined by the param- 
eter ta.b/o-Xi where ax — 5.4 nm is the exciton radius 
for CdSe.'' For the small QD's in Ref. 'E, the Coulomb 
interaction is small compared to the single-particle level 
spacing, such that the admixing of higher orbital levels 
to the ground state is small as well. (For details on ex- 
perimental parameters, see Sec. [3) This allows us to 
describe the coupled QD's by the Hamiltonian 



where 



Hn = 



E 

^A.B-a 



Coul 



(1) 



) (2) 



contains the single-particle levels of uncoupled QD's v = 
A, B. The operators cj^'"^ and c'^''^ annihilate an electron 
in the lowest level of the conduction band with spin 
quantum number Sz = crl/2 and the highest level in the 
valence band, E!^, with angular momentum jz — cr3/2, 
respectively, where a — ±. Here, we have adopted a 
simple model for the change in the band structure of 
CdSe due to the QD confinement. We assume a spherical 
QD shape and a splitting of the j = 3/2 valence band at 
the F point into the heavy hole (hh) and light hole (Ih) 
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subband with total angular momentum projection jz = 
±3/2 and jz = ±1/2, respectively, as obtained, e.g., from 
the Luttinger Hamiltonian with an additional anisotropy 
term for the crystal field of the hexagonal lattice. * The 
Ih subband will therefore be neglected in the following. 
The Coulomb interaction energy is 

U. 
2 



E 

'=A,B 



Hcon, = > : ^ KiK - 1) + KiK - 1) - 2n^<] 



where = Z^cr^-^ -c;^ , 



and = J^cr. 



(3) 

are the number operators for electrons in the conduc- 
tion band level and holes in the valence band level. 
Uu ~ /Ai:eeQri, is the characteristic charging energy of 
QD A and i?, respectively. Transfer of spin and charge 
between the QD's is accounted for by the transfer Hamil- 
tonian 



+ h.c.) , (4) 



where we assume that transfer of electrons through the 
TT-conjugated molecule conserves the electron spin both 
in the conduction and the valence band. 

The ansatz for the Hamiltonian in Eqs. is a 

model in which the biexciton shift, the exciton fine struc- 
ture, and the electrostatic coupling between the QD's 
have been neglected. We will justify this in Sec. |Vl be- 
low where we discuss our results for the experimental 
parameters of Ref. 5. Because the focus of this work is to 
calculate the FR angle that results from transfer of elec- 
trons between the QD's, we assume for simplicity that 
the symmetry axis of the QD's with hexagonal crystal 
structure is parallel to the direction of pump and probe 
laser pulses. The effect of a random QD orientation will 
be discussed in Sec. 

In the following, we analyze the results of Ref. phased 
on the Hamiltonian Eq. (^). This paper is organized as 
follows. In Sec. ^ we calculate the time-resolved FR 
signal for an electron wave function which is delocalized 
over QD's A and B. In Sec. lIIII we calculate the FR angle 
as a function of probe energy for an initial spin polariza- 
tion created by optical pumping. We take into account 
both electron transfer processes and the Coulomb inter- 
action and show that these terms give rise to an exchange 
splitting of the two-exciton eigenstates. In Sec. IIVI we 
perform the related analysis for a system with one spin- 
polarized excess electron in the QD's. In Sec. [3 we dis- 
cuss our results for the parameters of CdSe QD's coupled 
by benzene molecules;^ calculate the transfer matrix ele- 
ment and spin transfer probabilities. In Sec. I VII we draw 
our conclusions. 



II. TIME RESOLVED FARADAY ROTATION 
FOR COUPLED QUANTUM DOTS 

Before we calculate the FR angle for the general Hamil- 
tonian Eq. ^ in Sees. IIIII and IIVI below, we first con- 
sider time-resolved FR for a particularly simple case in 



which a single electron is in a coherent superposition 
of states in QD's A and B at time t = 0, \ip{0)) — 
(c^^. + Q!C^+)|0) /-\/l + a^. We further assume tc,v = 
and = E^ in Eq. ||TJ for t > 0. Here, |0) denotes 
the vacuum state in which the valence band in both 
QD's is filled and the conduction band states are empty. 
This simple scenario, although unrealistic because trans- 
fer matrix elements are assumed to vanish after the initial 
state I^/jq) has been prepared, will allow us to derive sim- 
ple analytical expressions for the FR angle even in pres- 
ence of a magnetic field. The simplifying assumptions 
tc.v — and E^ = E^ will be lifted in the microscopic 
discussion in Sees. IIIll and IIVI 

The different radii and of the CdSe QD's lead 
to different g-factors and different Larmor precession fre- 
quencies iOi, = giy ^J,B Bext / fij^'^'^'^^ where B^xt is an exter- 
nal magnetic field perpendicular to the spin quantization 
axis which is given by the symmetry axis of the CdSe 
QD's, and gi, are the electron g-factors for — A,B. At 
time t, 



\m) = 



1 



cos((jj_Bt/2)cf 1 — ism{iLjBt/2)cfl_ 



+acos{uJAt/2)cfl~iasm{ujAt/2)c^l |0). (5) 

This time evolution of the electron spin can be de- 
tected by FR because the FR angle 6p is determined by 
the population imbalance between the Sz = ±1/2 con- 
duction band states in this situation»iSiiii4 For probe 
pulse frequency E/h, Op is proportional to the difference 
of the real parts of the dielectric response functions €{E) 
for <7^ circularly polarized light With the spectral rep- 
resentation of the response functions, Op{E) is expressed 
in terms of the transition matrix elements between the 
state \ip{t)) with energy Eq and all intermediate states 
{ipi) which are virtually excited by the probe pulse, 

X (^\(i,,\p+\i,it))\^ ~\mp-\mf 

The polarization operators P± = dAcf]^c^^ + dBC^±c^± 
couple to the circularly polarized components of the 
probe pulse, di, are the dipole transition matrix elements 
for transition from the jz = ±3/2 valence band states to 
the Sz = ±1/2 conduction band states in QD's A and B. 
Eq = Ef and Ei are the energy eigenvalues of the initial 
state and the intermediate state respectively, and 
the level broadening F accounts for a finite lifetime of the 
orbital levels. The prefactor C cx L/ {hcno) is determined 
by the size L of the sample and the refraction index uq 
of bulk CdSe. 

Because we have assumed an initial state \tp{0)) with 
one electron, all intermediate states in Eq. ijHJl are 
energy eigenstates with two electrons and one hole. For 
tc,v = in Eq. QJ, these are of the form 1-0^) = 
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c^J^aCl^a^l'^^) '^ith cr,cr' = ± and v,v' = A,B. Pauli 
blocking prohibits the creation of an exciton with elec- 
tron spin (t1/2 if the conduction band level is already 
occupied by an electron with the same spin. The result- 
ing difference in transition matrix elements for P_|_ and 
P_ is proportional to the population imbalance of the 
Sz = ±1/2 levels. For a probe pulse at time i, from 
Eq. lO we obtain directly 
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E^y 



r2 



cos(cj^t) 



(7) 



where = E^ 



E'^ — Uu is the exciton energy for QD 



OpiE^t) shows coherent oscillations with frequencies 
LJA and lob caused by the electron spin precessing around 
the external magnetic field. In reality, these coherent os- 
cillations are exponentially damped with a spin dephas- 
ing rate Ts which is typically much smaller than the or- 
bital dephasing rate, T s ^ F. Taking into account spin 
dephasing, the Fourier transform of the time-resolved FR 
signal as a function of the probe pulse energy E and the 
Fourier frequency lj is 



OpiE^uj) 



CE 
l-f q;2 

E-El 



2j2 



E-E^ 



(8) 



6f{E,u;) shows characteristic features for E 
and a; ~ ijjj^. The two terms in Eq. (jS]) describe the 
dielectric response due to virtual creation of an exciton 
in QD A and B, respectively. For E^ < E < E^, they 
have different sign and may cancel. Figure ^a) shows 
a grayscale plot of \9p{E,llj)\ for the experimental val- 
ues E^ = 2.06 eV, E^ = 2.41 eV, F = 0.05 eV, and 
Fs/27r = 0.5 GHz, assuming d\/d% = 1 and = 0.2. 
For Fig.mb), F = 0.035 eV, and Fs/27r = 1.2 GHz, and 
= 0.4. One of the most characteristic features of the 
experimental data (Fig. 2C in Ref. 0) is that \9f{E,uj)\ 
vanishes and reappears as a function of probe pulse fre- 
quency E for Lo ~ lob- This feature is also present in 
the theoretical result and can be traced back to the su- 
perposition of two response functions in Eq. More 
specifically, for uj ~ ljb and E ~ E^ — F, the two terms 
in Eq. ((SJ have opposite sign and cancel for sufficiently 
large a. 

Above, we have assumed that the electron delocalized 
over both QD's at t = retains spatial coherence. For 
rapid decoherence of the orbital part of the wave func- 
tion, the initial state is described by the density matrix 



P 



St 

+ 



a c, 
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FIG. 1: (a) Grayscale plot of the FR angle \9f\ [given in 
Eq. (|SJ] in arbitrary units as function of probe pulse frequency 
E/h and frequency / = uj/2n. We have chosen the param- 
eters^ = 2.41 eV, S| = 2.06 eV, lua/2-k = 23.6 GHz, 
ujb/2-k = 20.6 GHz, V = 0.05 eV, rs/27r = 0.5 GHz, and 
= 0.2. (b) For F = 0.035 eV, rs/27r = 1.2 GHz, and 
= 0.4, pronounced features caused by the interplay of the 
two terms in Eq. Q become more clearly visible. In partic- 
ular, the FR signal at uj ^ uib vanishes and reappears as a 
function of probe pulse frequency E. 



FR signal in this case is the incoherent superposition of 
the FR signals for QD A and B, and is identical to the 
results in Eqs. Q and ©. Hence, a FR signal as shown 
in Fig. n does not allow one to distinguish coherent from 
incoherent spatial superpositions. 



III. OPTICAL SPIN INJECTION 

In the preceding section, 0f{E) was calculated for the 
simple case of a single electron delocalized over the cou- 
pled QD's. So far, we have also neglected that all inter- 
mediate states IV'i) in Eq. © that are virtually excited by 
the probe pulse will be modified by finite transfer ener- 
gies tc,ii- We next turn to a microscopic analysis in which 
we take into account tc^v 7^ also for the intermediate 
states. 
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In Ref. 0, the initial state prepared by optical pump- 
ing is a one-exciton state. Similar to the analysis in 
Sec. m above, the FR angle as a function of probe 
energy is proportional to the difference of the dielec- 
tric response functions for circularly polarized light 
[Eq. jnj]. In order to evaluate this expression, both the 
initial one-exciton state and all intermediate two-exciton 
states which are virtually excited by the probe pulse must 
be calculated for the coupled QD's. In this Section, we 
first calculate the one-exciton energy eigenstate of the 
coupled QD's prepared by the pump pulse and subse- 
quently identify all two-exciton eigenstates \ipi) which 
are virtually excited by the probe pulse. Our analy- 
sis is based on perturbation theory in the transfer en- 
ergies and is valid if \tc,v\ is the smallest energy scale, 
\tc^y\ < 5Ec, \5Ey\, Ua, Ub, \SEc^y±UA,B\- Here, we have 
defined the energy differences SE^ = Ef - E^ > Q and 
5Ey = E^ — E^ < between the conduction and valence 
band levels of QD's A and B. 

In Ref. 0, an initial spin polarization was created by 
optical pumping. For tc,v = 0, the states cj^'l'^c^ g.|0) are 
one-exciton eigenstates with energy eigenvalues 



E 



X 



Er. - Ef, 



(9) 



which are prepared by absorption of a —cr circularly po- 
larized pump pulse. To first order in the transfer energies 
tc^v, the energy eigenstates are 



\Xa,.) = c^^<jO)-|-(^— 



Ua 



SE, + Ua 



c..a |0), (10a) 



tc 



5Ec + Ub 



ty 



SE, - Ub 



cfX^^yO), (10b) 



with eigenenergies 



E 



A(0) 
X 



ii 



ti 



SEc - Ua SE., + Ua 



Ei = 



ti 



L,,. 



5E, + Ub SE, ~ Ub 



(11a) 
(lib) 



We now turn to the calculation of the FR angle, as- 
suming that the pump pulse has prepared an initial state 
lip) = \Xb,+)- The evaluation of the dielectric response 
function will require us to calculate all two-exciton states 
that are virtually excited by the probe pulse. Interesting 
features in the FR signal effected by spin transfer are of 
order t'^ ^. In order to keep the following expressions sim- 
ple, we assume that spin is transferred between the con- 
duction band states and set ty — 0. Then, only the seven 

states \A+B+), |To), \S), \B+B-), |Tb), \S), and jS^^) 
listed below and in Appendix ^ have finite matrix ele- 
ments up to 0{tl) with P±\Xb,+). For SEy + Ua 0, 
only the eigenenergies of \A+B+), |To), and \S) are close 
to the excitation energy of a probe pulse with frequency 
E/h ~ E^/h. Hence, these states dominate the spectral 
representation in Eq. 

The polarization operator P_|_ induces transitions from 
the initial state \Xb,+) to 



\A+B+) = cilcfU^^c: 



|0) 



with energy eigenvalue 



Ea^b^ = E 



A{Q) 
X 



E 



B{a) 

X 



(12) 



(13) 



The notation indicates that two electrons with the same 
spin Sz = 1/2 occupy the conduction band states in 
QD's A and B, respectively, and form a spin triplet state. 
\A^B^) is an exact eigenstate of the Hamiltonian even 
for ic 7^ because transfer of the conduction band elec- 
trons is blocked by Pauli's exclusion principle. The ma- 
trix element {A+B+\P+\X Bj^) is the only finite matrix 

element of the operator P+ . 

Finite matrix elements for P_ come from the states in 
which the electrons in the conduction band level form a 
spin triplet and singlet, respectively, 

l^o) = ^ (e^f I + ^clcT-) <-C^+|0), (14a) 

\s) « i= (cf_af^ _ a^^af ^) c^cf,+|o) (i4b) 

+%/2 



tr 



tr 



SE, + Ub 5Ec - Ua 

xc;f^_c^+|0). 



As expected, the eigenenergies are shifted due to the 
delocalization of electrons and holes over the coupled 
QD's. The exciton states in Eq. (|l()|l are the only one- 
exciton states which can be prepared by the absorption 
of a photon with circular polarization —a if the photon 
is incident along the hexagonal axis of the CdSe crystal 
structure. However, a photon with energy E ~ E^ no 
longer creates an exciton only in QD B, but an exciton in 
which electron and hole are delocalized over the coupled 
QD system. This delocalization of the quantum mechan- 
ical wave function is consistent with the short time scale 
for spin transfer observed experimentally*^ 



and the holes with jz = —3/2 and jz = +3/2 are local- 
ized in QD's A and B, respectively. Note that the pro- 
jection of the total conduction band spin onto the spin 
quantization axis vanishes for the triplet state iTo)- The 
normalization constant for \S) is defined by {S\S) = 1. 
The eigenenergies 
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FIG. 2; Schematic representation of the spin configurations 
(in the electron picture) for states (a) \A+B+) and (b) IS), 
|To) to leading order in tc- The dashed hues represent the 
conduction and valence band edge in bulk CdSe. 



show an energy offset which is caused by the inter-dot ex- 
change coupling. "''^'^ The energies of \A+B+) and \Tq) are 
not shifted by electron transfer because of Pauli block- 
ing and destructive interference of transfer paths, respec- 
tively. 
The state 



|B+B_) cx 



tr 



SE, - Ua 



{5E,-Ua){25E,-Ua + Ub) 



-At -St _ --41 -St 



ct-cl^\0) (16) 



with 
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B(0) 
X 



Ei 



E. 



5E, - Ua 



(17) 



is offset in energy from E^'"'^^ 



-'X 



-|- £'^*-°-' even to zeroth or- 



der in tc and does not contribute significantly to Of{E) 
for E ~ E^. The three states in Eqs. ^ and ^ 
provide the dominant terms in the spectral representa- 
tion for dp in Eq. ©. In particular, they exhaust the 
sum rule | (V-d^Lc^,. |Xb,+) p = 1 up to Oitl). In 

Fig. 13 the spin configurations for \A^B+), \S), and \To) 
are shown schematically. 

From Eqs. l|10b(l - (fT7|l . the FR angle 0f is readily eval- 
uated. We denote the electron transfer probability from 
QD 1/ to QD ly' by p^^^i . We obtain 



PA^B 



Pb^a = 



tc 



SEc-UaJ ' 
t ^ ' 
5Ec + Ub 



(18a) 
(18b) 



For the transition matrix elements of the dipole operators 
in Eq. (j^ , we obtain in terms of the transfer probabilities 



\{A+B+\P+\Xb.+) 
\{To\P-\Xb,+)\' = 

\(S\P^\Xb,+)\' = 



{l-pB^A)d\, (19a) 
(19b) 



l_PB^^, 



1 + PB^A - "^.PA^B j2 nn \ 
n ^A, (19c) 



{B+B^\P.\Xb,+)\^ ^PA^Bd\. 



(19d) 



Because of the exchange splitting Etq — Es between 
conduction band triplet and singlet states, finite transfer 
probabilities pa~*b and pb-^a lead to pronounced fea- 
tures in the FR angle as a function of the probe pulse 
frequency E/h. For probe energies EtqB — Et„ — E^ < 
E < Es B = Es — E^ , the FR signal varies strongly with 
energy and is given by 



Of{E) = 



CE(f, 



(1 ~pb^a) 



E-E. 



TqB 



(1 +pb^a - '^pa^b) 



(E-ET,B)^ + r2 
E — EsB 



{E-EsBV+r\ 
FoY \E - Esb\>\Et„ - Es\, Eq. 12011 simplifies to 



(20) 



CEdl 
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A(0) 
X 



^X 



{PA- 



>S - Pb^a) 



Y2 



(21) 

This result is surprising because the FR angle is not only 
determined by the probability pb^a that the electron 
created by the pump pulse has been transferred to QD A. 
Rather, even the sign of the FR angle depends on the pa- 
rameters 5Ec (and 5E^ if transfer between valence band 
states is included) and Ua,b- Of > for \5Ec — Ua\ > 
\5Ec -h f/s I, and < for V^^c - Ua\ < \5Ec + Ub\- Al- 
though counterintuitive at first sight, this can be readily 
understood from the one- and two-exciton eigenstates. 
The matrix element for the virtual creation of an exci- 
ton with Sz = 1/2, jz = 3/2 in QD A is reduced by the 
probability that the conduction band electron cre- 

ated by the pump pulse in B has been transferred to A. 
In this case, it blocks the creation of a second exciton 
with the same spin. The transition matrix element for 
the creation of an exciton with Sz = —1/2, jz = —3/2 is 
reduced by the probability pa—^b that the electron with 
spin Sz = —1/2 in the conduction band state of QD A is 
transferred to QD B. This transfer process is not prohib- 
ited by Pauli blocking and leads to the virtual occupation 
of jS+B-) which is energetically far off resonance. The 
interplay of both processes results in Eq. (|21|l . 

Our derivation of Eq. (|21|l was based on the assump- 
tion that tc is the smallest energy scale in the system. As 
will be discussed in Sec. below, for the experimental 
parameters in Ref.^^ SEy^UA — 0. For ty — 0, this does 
not lead to divergencies in the perturbative expansion in 
tc- However, these special parameters require that two 
additional two-exciton states are taken into account for 
the calculation of 6f{E) because they are nearly degen- 
erate with \A+B+), \S), and \To) (see Fig.OJ: The states 



6 



pA(0) , pB (02 

l_y I I— Y 



2EB(o) __ 



IS) 



|A,B,>, |To) 

it;) 

IS) 

|B,B.) 
I^B ) 



:}o(tc 



■}0(5EJ 



5E, + U, 



0(5E, 



FIG. 3: Energy level scheme of all two-exciton eigenstates 
discussed in the text. The eigenenergies fall into three groups 
which are split by terms of order 0{tc) or 0(5Ev). For the 
QD's used in Ref.|3, SE^ + Ua — 0, and the five states \A+B+), 
I To), 15"), I To), and 15") are nearly degenerate. 



\S) and |To) defined in Eq. lAip have finite overlap ma- 
trix elements with P-\Xb,+), 



\{fo\P-\XB, + ) 

\{S\P^\Xb,+) 



2 ' 



(22a) 
(22b) 



The spin configurations for the states \S) and \To) are 
shown schematically in Fig. QIa). Note that both holes 
occupy the valence band states of QD B. The accidental 
degeneracy of \S) and \To) with \S) and |To) arises be- 
cause, for the parameters of Ref.0, the decrease in orbital 
energy 6Ey is comparable to the increase in Coulomb en- 
ergy Ua- Transitions between an initial state \Xb,+) and 
|S'), \Tq), are two-step processes. A (t+ polarized probe 
photon creates an exciton with = — 1/2 and jz — —3/2 
in B, and one of the conduction band electrons in B 
is subsequently transferred to A. These processes are 
shown schematically in Fig.^fb). 

Taking into account all two-exciton states with ener- 
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The energy differences Ej, g 

Eg — E^ are given by the eigenenergies in Eq. (|A2|) . For 
\E-Et,b\AE-Esb\ « \E-Ef^gl\E-EsglEq. ^ 
simphfies to Eq. if^ 
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'0 /' 



s = 



+ 



t- 



(b) 







■■■■■^ 



FIG. 4; (a) Schematic representation of the spin configura- 
tions for the states \S), |To) to leading order in t^. (b) Transi- 
tions between an initial state \Xb,+} and \S), |To) are effected 
by the absorption of a a'^ polarized probe photon and subse- 
quent tunneling of one conduction band electron. 



Above, we have only considered tc and ty — 0, 
i.e., a scenario in which electrons in the valence band re- 
main localized in the QD's while electrons in conduction 
band states can be transferred. The case ty and 
tc — can be mapped onto the problem discussed above 
by mapping electrons onto holes, i.e., by interchanging c 
and V in above expressions. In particular, Eqs. (|2U|) and 
(|21|l remain valid if the transfer probabilities for electrons 
are replaced by the corresponding values for holes, e.g., 
PA^B = [tv/{SEy + Ua)]^, and the energy eigenvalues 
are calculated for transfer in the valence rather than the 
conduction band. 

In the limit of small QD's with similar sizes, Ua.b ^ 
tc,v S> SEc, \6Ey\, configurations in which electrons and 
holes occupy different QD's are strongly suppressed. If 

tc,v/UA,B ^ but tcty/UA.B (^eP°^ - Ex'^^^'j remains 
finite, a joint transfer of electron and hole via a virtual 
intermediate state is possible. Evidence for this coherent 
delocalization of an exciton has been reported for QD's 
of similar sizesjiS^i^ In contrast, tunneling of excitons be- 
tween a pair of QD's with different sizes is incoherent if 
the orbital decoherence rate is comparable to the exciton 
tunnehng ratei2Si2ii22iS^ 



IV. DOPING OF COUPLED QUANTUM DOTS 

In the last Section, we have analyzed the FR angle for 
an initial spin population created by optical pumping, 
the method used in Ref. 0. We now calculate the FR 
angle 9f{E) for the case that the initial spin density is 
carried by an excess electron rather than an exciton. Spin 
injection could be achieved, e.g., by doping one CdSe 
QD with a single donor atom. For a chemical potential 
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< IJL<E^, E^+Ub, the conduction band level of QD 
B is filled with one electron while QD A remains empty. 
The excess electron can be spin polarized by cooling in 
presence of a magnetic field. Again, we set i„ — to keep 
our results transparent. 

The transfer matrix element for the conduction band 
level leads to the delocalization of the excess electron in 
QD B, 



eB,cr) = 



tc 



-1/2 



tc 
SEr_ 



ctl]\0) (24) 



with eigenenergy = Ef — tl/SEc- Note that the 
energy shift is different from the one found for the exciton 
because there is no Coulomb attraction between electron 
and hole in the present case. 

We calculate the FR angle for an initial state \eB.+) 
and probe energy E ~ E-^. Similar to the analysis in 
Sec. mil three intermediate states dominate the spectral 
representation for 6f{E). These states are the following. 



eigenstates with two conduction band electrons and one 
hole listed in Appendix ^ are energetically far off reso- 
nance and can be neglected. 

The transition matrix elements of the polarization op- 
erators P± between |eB^+) and the states Eqs. (|25() and 
(|27|l are readily evaluated. The probabilities for electron 
transfer between the QD's are now given by 



PB-.A 

Pa^b 



_tc_ 

6Er 



tn 



5Ec -Ua-Ub 



Then, 



\{T,\P-\eB,+)\' =^-^^d\, 



(29a) 
(29b) 



\{A+Bl\P+\eB,+)\' = {l-PB^A)d\. (30a) 

(30b) 



\{S-\P4eB,+)? = i±P^^4-^^4.(30c) 



\A+B-)=c^yc'!Xci+\Q) 
with energy eigenvalue 



Ea+b- - Ex +Ec 



(25) 



(26) 



is populated by creation of an exciton with conduction 
and valence band spins — 1/2 and jz ~ 3/2, respec- 
tivelyi24 Virtual creation of an exciton with Sz — —1/2 
and jz — —3/2 leads to transitions to the spin triplet and 
singlet states 



To) = ;^(ee+ee)fiMo) 

1 / 



(27a) 
(27b) 



'^^Ui?c''^'+''^'- SEc- Ua-Ub 

xc;^.-|0), 



where the normalization constant for \S ) is determined 
by (S'~|S'^) = 1. The eigenenergies 



Err^- 
^0 

E^- 



eP"^ + e^, 

E^^"'' + E^ 



(28a) 
(28b) 



+2ti 



1 1 

SEc ~Ua-Ub^SEc 



Inserting these matrix elements into the spectral repre- 
sentation of 9p{E), Eq. lO, we find for the FR angle 



OAE) = ^ 



,^-Pb^a> 



E~ E. 



T„B 



{^+Pb^a-^P 



(E-E-^^r + r^ 

E - 



A^BJ 



for probe energies E ~ E^'"^^ , in close analogy to Eq. H20() 
for optical spin injection. The energy differences are de- 
fined by E^ g 



{E-EsBr + r^ 



(31) 



Erp- 



E^ and E 



SB 



Es- - E". Be- 



cause of the exchange splitting between iTg") and 15*"), 
6p(E) will in general exhibit several peaks and lack point 
inversion symmetry. The functional dependence on probe 
energy is determined by the transfer probabilities and the 
energy differences E^^^ and Egg. For a more detailed 
analysis which takes into account all finite transition ma- 
trix elements up to 0{tl), see Appendix IBl 

Experiments on doped QD's could provide valuable in- 
formation supplementing the experimental data obtained 
for optical pumping. The main advantage over optical 
spin injection is that spin decoherence times are expected 
to be substantially longer because they are not limited 
by electron-hole recombination. Even more importantly, 
FR measurements on doped coupled QD's can clarify 
whether spin transfer occurs predominantly between the 
conduction or valence band levels because, for tr = and 



ty ^ 0, 9f{E) ~ for probe energies E ~ E 



A{Q) 

X ■ 



are split by the exchange coupling of the conduction band 
levels. Further, there are several states with energies dif- 
fering from E^^^^ + Ef (see Appendix . For probe 

pulse energies E ~ £'^'"^ and \5Ey 4- [/^ - ^ F, 
9p{E) is dominated by virtual excitations into the states 
\A+B^), [Tg ), and \S^). In this case, all other energy 



COMPARISON WITH EXPERIMENT 



In order to compare the results of Sec. IIIII with ex- 
perimental data from Ref. |^ we first provide numeri- 
cal values for 6Ec, SE^, Ua, and Ub- The energy level 
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spectrum of CdSe QD's is well established both experi- 
mentally and theoreticallyiSSi2& The absorption energies 

E^^°^ = 2.41 eV and Ep°'^ = 2.06 eV in Ref.Hare consis- 
tent with VA — 2.0 nm and rs — 3.5 nm, and we will use 
these radii for the following calculations. From Ref. |2^ 
6Ec ~ 0.30 eV and 5E^ ~ -0.10 eV. 

From the bulk values for the static dielectric constant, 
e = 9.7, and the band masses in the conduction and 
valence band, mc/rrif^ — 0.12 and my/mi, — 0.45, one 
obtains the exciton radius 5.4 nmiSLSSiSS The exciton 
radius is larger than r^,_B, and electrons and holes are 
strongly confined in the QD's as assumed in Eq. The 
characteristic energy scale of the Coulomb interaction is 
Ui, ~ e^/47reeory. For the given values of rA and rs, 
Ua = 0.07 eV and Ub = 0.04 eV. 

The Hamiltonian Eq. does not take into account 
biexciton shifts, the exciton fine structure, and inter-dot 
Coulomb interactions. For CdSe QD's with radii 1.5- 
4 nm, the biexciton shift is of order 0.01-0.02 eV fRef. l30|) 
and the characteristic energy splitting between bright 
and dark excitons is smaller than 0.01 eVj2i The char- 
acteristic energy scale for inter-dot Coulomb interactions 
is Uab — e^/47reoe(r^ + rs) < 0.03 eV. However, it is 
relevant only if neither of the two QD's is electrically neu- 
tral. The most important effect of the inter-dot Coulomb 
interaction is to lower the energy eigenvalues of |To) and 
IS") [Eq. HA1|) ] by Uab- All these energy scales are small 
compared to the level broadening F and can safely be 
neglected. 

In the following, we assume that only electrons in con- 
duction band levels are transferred between the QD's 
while valence band electrons remain localized. As dis- 
cussed in Sec. IIVI this assumption can be tested by ex- 
periments on doped QD's. Mediated by electron trans- 
fer through the molecular bridge, the lowest conduction 
band level in QD B hybridizes with the lowest conduc- 
tion band level in QD A. Comparing the observed energy 
shift E^ ~ E^^°^ = -0.02 eV with Eq. (EJ, we find 

tc = ^ (-^x'"^ - {SEc + Ub) = 0.082 eV. (32) 

Our theory predicts that the exciton absorption peak 
for QD A is shifted to larger energies for the cou- 
pled QD's, in contrast to the experimental result E^ — 
E^^^^ < 0. The most likely explanation for this is 
that the lowest conduction band level in QD A hy- 
bridizes also with higher excited levels in QD B which 
are nearly degenerate with E^fi^ In order to account 
for quantitative changes effected by this hybridization, 
the energy i?^^"^ must be replaced by the true value 
of the hybridized state in all expressions for the two- 
exciton eigenenergies. This value can be obtained from 

■E^x °^ + tl/{SEc + Ua) ^ 2.36 eV, where the latter is the 
experimental value for the exciton absorption edge of QD 
A in the coupled QD's. Hence, E^^^ 2.33 eV. 

From these parameters, we calculate for the transfer 
probabilities between the lowest conduction band states 



PA^B = 0.13 and pb^a = 0.06. The energy differ- 
ences between the two-exciton states and the initial state 
are Et„b = 2.35 eV, Esb = 2.37 eV, Ef^g = 2.32 eV, 
Egg — 2.31 cV. The oscillator strength for exciton cre- 
ation, proportional to d\ is independent of the QD 
size in the strong confinement regime and proportional 
to the QD volume for weak confinement. Because both 
QD's are close to the strong confinement limit, we assume 
a weak scaling c?^/c?^ = 2 for the following Figures. 

In Fig. EJa), we show the FR angle calculated from 
Eq. (|21|l as a function of probe energy for different values 
of r, F = 0.05 eV (solid), 0.02 eV (dashed), and 0.08 eV 
(dotted). We note that even quahtative features depend 
strongly on the microscopic parameters such as F. For 
small F, additional peaks emerge because the contribu- 
tions from the individual two-exciton states can be re- 
solved. 

In spite of the dependence on microscopic parameters, 
some pronounced features in 6p(E) are generally present: 
(i) 9f{E) does not exhibit point- inversion symmetry, in 
stark contrast to the FR angle expected from virtual 
transitions to a single state, (ii) Op has in general more 
than two maxima or minima. The positions and heights 
of the extrema are determined by the interplay of the 
transfer probabilities pa^b and pb^a, and the energy 
splittings between the different two-exciton states. Ex- 
periments have demonstrated the strong dependence of 
the FR angle on the probe energy E, including a fine 
structure of the resonance. '^■^ 

In Fig. [SJb), we compare the calculated FR signal for 
coupled QD's A and B with the corresponding result for 
uncoupled QD's A pumped at resonance. For a probe 
energy E ~ 2.42 eV, the FR signal for coupled QD's A 
and B is significantly smaller than the FR signal of the 
AA system, consistent with experimental observations?^ 

Note that in Fig. |S1 we show the FR angle in arbi- 
trary units because its absolute numerical value depends 
sensitively on unknown experimental parameters such as 
the packing density of QD's and the number of QD's A 
which are coupled to at least one QD B. For a spin 
transfer probability pb—>a = 10%, assuming close pack- 
ing of the QD's and that every QD A is coupled to one 
QD B, we estimate the maximum value of the Faraday 
rotation angle to be 0.02 rad for an ABAABA-atTuctme 
as investigated in Ref. 0. 

So far, we have assumed that the symmetry axis of the 
CdSe QD's with hexagonal crystal structure is parallel to 
the propagation direction of pump and probe laser pulses. 
However, in experiment the QD's are randomly oriented. 
We discuss next how the random orientation changes our 
results. The propagation direction of pump and probe 
laser pulse is z, the polarization vector of the probe pulse 
X, and the symmetry axes of QD's A and B are denoted 
by ca and cb, respectively. We define the azimuthal 
angles (pA = ^(x,ca) and (j)B = Z(x,cb), and the angle 
enclosed by the two symmetry axes (pAB = ^{caiCb) 
[see Fig. EJa)]. The conduction band spin eigenstates 
with quantization axis ca,b are denoted by | Ta,b) and 
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(a) 



/ \ 
/ \ 






V ^ — 



2.25 2.3 2.35 2.4 2.45 
E[eV] 

(b) 




^ ^ ^ _ ^ 

2.25 2.3 2.35 2.4 2.45 
E[eV] 

FIG. 5: (a) Plot of the FR angle as a function of probe pulse 
frequency calculated from Eq. H21|l for different level broad- 
enings F = 0.05 eV (solid), 0.02 eV (dashed), and 0.08 eV 
(dotted). All other parameters are as described in the text. 
For small F, 9f{E) clearly shows the individual contributions 
from the various two-exciton states, (b) Comparison of the 
FR angle for coupled QD's for F = 0.05 eV (solid) with the 
calculated signal for a AA structure (dashed). 



I Ia^b)- 

For arbitrary angle Z(z, c^), the probability for the 
circularly polarized pump pulse to create a net spin po- 
larization in the conduction band level decreases from its 
maximum value at Z(z,cb) = to zero at Z(z, c^) = 
tt/2. For Z{z,cb) < tt/2, the majority of conduction 
band electrons is in spin state | ts), with the quanti- 
zation axis defined by c^. On transfer to QD A, the 
conduction band electron retains its spin state because 



states with Sz = ±1/2 are degenerate in both QD's and 
tc is spin- independent. The characteristic level spacing 
of valence band states is large compared to the crystal 
field splitting in bulk CdSe, which allows us to treat the 
latter as a small perturbation, following Ref. [s^- 

In the following, we calculate the FR angle for a ran- 
dom orientation of QD's assuming that the pump pulse 
has created a conduction band electron with spin | ts)- 
The random orientation of QD's affects the FR of the 
probe pulse in two ways. Firstly, the matrix elements 
for transitions from the jz = ±3/2 valence band lev- 
els to the Sz = ±1/2 conduction band levels in QD A 
(B) decrease by sini/)^ (sin(/)s) compared to the oriented 
sample. More importantly, also the relative orientation 
of ca and cb modifies the FR angle. For illustration, 
consider two QD's with tc — 0, and a conduction band 
electron in spin state | ts) in -B. The cr^ circularly 
polarized component of the probe pulse with E ~ i?^ 
excites a virtual exciton in A, with a conduction band 
electron in spin state | 1a)- Note that the spin direction 
is defined by ca, the symmetry axis of A. Expanding 
I 1a) = cos((/)ab/2)| 1b) ± isin(0^B/2)| Ib) in terms of 
the eigenstates along quantization axis c^, the product 
state of the two excitons contains terms in which the two 
conduction band spins are antiparallel and have a finite 
overlap with the spin singlet state. This is in stark con- 
trast to the oriented sample, where the two conduction 
band electrons would always form a triplet. 

The analogous analysis for coupled QD's must take 
into account both the reduced transition matrix elements 
for the probe pulse and the relative orientation of QD's 
A and B. Because virtual transitions to |To) and 
involve excitation of QD B which was populated by the 
pump pulse, the matrix elements in Eq. (|22l) are reduced 
by a factor | sin 0s | which is independent of the relative 
orientation of ca and c^. In contrast, virtual transitions 
in QD A probe the spin polarization relative to the quan- 
tization axis Ca after an electron with spin pointing along 
Cb has been transferred, and the transition matrix ele- 
ments depend also on (pAB [Fig-E^b)]. For the FR angle, 
we find 



9f{E) = { d\ cos 0AB sin^ (j)A 



{1-Pb^a) 
d| sin^ (j)B Pb 



E - EtoB ,^ o \ E - EsB 

{E-Er.B)' + V^ - ('+P^-^ - iE~EsBr + r\ 



E-E. 



TnB 



E-E 



SB 



{E-E, 



ToB) 



F^ [E-Esb)^ + T^ 



(33) 



The dependence on the relative orientation of the two 
QD's, is readily understood. For 0ab — 7r/2, the 

first and second term in the expression for Op {E) vanish 
because the conduction band spin created in QD B is 
perpendicular to the spin quantization axis in QD A. A 



laser pulse probing QD A does not show any FR because 
the net spin along ca vanishes [Fig. El^b)] . 

In experiment, ca and cb are randomly distributed 
over the unit sphere. Performing this average in Eq. 
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QD A 




FIG. 6: (a) The hexagonal symmetry axes of QD's A and B 
are in general oriented randomly relative to the direction of 
the laser pump and probe pulse. Because of the interband 
selection rules, a a" circularly polarized laser pulse generates 
a spin polarization along the symmetry axis of the respective 
QD. (b) A conduction band electron created in QD B retains 
its spin direction on transfer to QD A. FR in QD A probes 
the projection of this spin onto the symmetry axis ca, which 
gives rise to a factor cos (pAB for the first and second term in 
Eq. CT . 



we find for the FR angle 



CE_ r_3_ 

2 \l6^ 



(1 - PB^A 

E-E, 



E — EtoB 



{E-Et,b? + T^ 
E — EsB 



{E-EsB? + T^ 



TqB 



(£;-i?^„)2 + r2 



+ 



^SB 



{E-EsB? + r' 



(34) 



Note that the spectral weight of the last term increases 
compared to the oriented sample. 



VI. CONCLUSION 

We have calculated the Faraday rotation angle for cou- 
pled QD's as a function of the probe pulse frequency. 
We have considered an initial spin polarization in neu- 
tral QD's (created by optical pumping) and of one excess 
electron in the two coupled QD's. Our results lead us to 
the following conclusions. 

(i) The Faraday rotation angle shows a nontrivial func- 
tional dependence on the probe energy, the details of 
which depend on the spin exchange energy and spin 
transfer probabilities [see Eq. (|21|) and Fig. I^^a)]. Most 
notably, because several two-exciton states are separated 
in energy by a small spin exchange coupling, 6p{E) is not 



invariant under point inversion symmetry. Measurement 
of 9f{E) as a function of probe energy would allow one 
to identify the contributions of the various two-exciton 
states that are virtually excited by the probe pulse. 

(ii) Experiments on doped QD's would allow one to 
determine whether spin transfer is mediated by transfer 
in the conduction or valence band states. In particu- 
lar, from a vanishing Faraday rotation angle for probe 
pulse energies close to the resonance of QD A one could 
exclude that an excess electron injected into QD B has 
been transferred to A. In contrast, for optical spin injec- 
tion, spin could be transferred both between conduction 
and valence band states. 

(iii) In general, measurement of the Faraday rotation 
signal at a given probe frequency does not provide enough 
information to determine spin transfer probabilities be- 
tween the QD's. However, from the experimentally ob- 
served energy shifts, we calculate a characteristic energy 
scale tc = 0.08 eV for spin transfer in the conduction 
band. Based on the transfer Hamiltonian ansatz, this 
implies a probability of 6% for electron spin to be trans- 
ferred from QD B to QD A, and of 13% for the opposite 
direction. 

The purpose of this work was to establish the connec- 
tion between spin transfer and the Faraday rotation sig- 
nal observed in experiment. Our analysis was based on a 
transfer Hamiltonian ansatz. Some of the most interest- 
ing results of Ref. remain to be explored theoretically. 
Most notably, the transfer Hamiltonian ansatz is based 
on the assumption that electrons are transferred between 
the QD's via the bridging benzene molecule. Microscopic 
work will have to clarify why conjugated molecules pro- 
vide efficient transfer paths between QD's. 

The results obtained here can provide important guid- 
ance also for the identification of microscopic transfer 
mechanisms. The increase of the Faraday rotation sig- 
nal at a fixed probe frequency has been interpreted as 
increase of the spin transfer efficiency for higher tem- 
peraturesn^ According to our results, an increase in the 
transfer matrix element tc also leads to a shift of the ex- 
citon edge in absorption spectra toward lower energies. 
If the exciton absorption edge does not change with in- 
creasing temperature, the increased Faraday rotation sig- 
nal is more likely effected, e.g., by additional incoherent 
transfer paths than by an increase of the transfer matrix 
element. 
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APPENDIX A: TWO-EXCITON EIGENSTATES 
OF COUPLED QUANTUM DOTS 

In order to evaluate the FR angle 0f{E) from Eq. 
for arbitrary probe energies E, all two-exciton interme- 
diate states \ipi) with finite transition matrix elements 
{^pi\P±\XB^+) must be calculated. States with energies 

Ei ~ E+E^^^^ lead to the dominant contributions in the 
expression for the FR angle, Eq. (jSJ. The states \ A+B^), 
\S), and I To) defined in Eqs. ifT^ and ((Tl|l have energy 



eigenvalues Ei with \E^^°'' +E 



pB{0) 
X 



E,\<0[tl/{SE, 



UA),t2/iSEc + Ua)], and are the most important inter- 
mediate states for probe pulse energies Ec:i E 



A{0) 
X 



How- 



ever, for the experimental values of Ref. ^ 5Ey + Ua 
is small and two additional two-exciton states must be 
taken into account. 
The states 



V2 



^l£^X + ^i^^^.+^.-\^, (Ala) 



n/2 V 
-V2 



tr 



SEc + Ua + 2Ub 
"SE, 



(Alb) 



differ from the corresponding states in Eq. (|14|1 in that 
both holes are localized in QD B. The normalization 
constant for \S) is fixed by (5^15) = 1. The eigenenergies 



with 



ti 



6E, + Ub 



(A4) 



is energetically separated from Et^ and Es by i?^^"^ — 



E 



B{0) 
X 



In Table HI we summarize all two-exciton eigenstates 
which contribute to the spectral representation of 9p{E) 
up to order t^. We also list the formal expressions for 
their eigenenergies to leading order 0{t^^) and give the 
numerical values, taking into account terms up to 0(t^) 



1^0 




E, - E§ [eV] 


1^+5+) 




2.35 


|To> 


^MO) _|_ ^S(O) 


2.35 


\S) 


£.-4(0) _|_ ^B{0) 


2.37 


\B+B-) 




2.06 


\fo) 


E^-E^ + E"/'^ 


2.32 


\S) 


E-^ - E^ + S^'"' 


2.33 


\B+B-) 




2.04 



TABLE I: Two-exciton eigenstates which contribute to 
the FR angle up to second order in tc- We also list the cor- 
responding eigenenergies to 0{t^) and evaluate them for the 
parameters discussed in Sec. El As noted in the main text, 
the degeneracy of \To), \S) with |To), IS) is a consequence of 
SEv + Ua — for the QD's used in experiment. 



Erf 



E 



X ^ ^x 



Ei ^ E 



A(0) 



3(0) 
X 

^eT^ 

+2t 



6E, 

SE„ 
1 



Ua, 
Ua 



1 



SEc + Ub SEc 



are shifted relative to Etq and Es by SE^ 
The state 



+ 



tr 



2tiCc]+C, 



SEc 

At 



Ub 



(A2a) 
(A2b) 



2C/bH 
+ Ua. 



Ua 



{SEc + Ub){2SEc + Ua + 5Ub) 



|0) (A3) 
I 



APPENDIX B: EIGENSTATES OF DOPED 
COUPLED QUANTUM DOTS 



Here, we calculate eigenstates and energy eigenvalues 
for states with two electrons and one hole in the coupled 
QD's. These are the intermediate states j-i/'i) in Eq. © 
which have finite overlap matrix elements with P±\eB.+) 
and determine the FR angle for coupled QD's doped with 
a single excess electron. 

In addition to the states \A+B^), \S^), and \T^) de- 
fined in Eqs. (|25|l and H27|) . five states have contributions 
of order to the FR angle. These are 
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\A:B. 



tc 



~Ua-Ub 



.,-10), 



1^0-) - 
\~S-) (X 
isTsl ) cx 



1 



%/2 



-At -St , 



1 /-^ 



?.A\ ^St _ -At -St 



V2 



tr 



-St -St 



gAt -St _ -At gSt 



5Ec + Ua + Ub 



^A\ -At 



-St -St 



r 



(Bla) 

(Bib) 
(Blc) 

(Bid) 

(Ble) 



with the proportionahty constants chosen to ensure nor- 
mahzation. The corresponding energy eigenvalues are 



^B+BZ 


= 2i;f + 


Ub - E^ 


(B2a) 







''c 


(B2b) 




5E, -Ua~Ub' 


E- 

A+B+ 






(B2c) 


^0 






(B2d) 




= _ 


VEt 


(B2e) 




+2tl 


i' - 






\SE, 5E, + Ua^ 


hC/s) ' 


E 

B+B- 


= E^^"^ - 




(B2f) 



From Eq. (|B1|I , we obtain the transition matrix elements 
in terms of the transfer probabilities defined in Eq. (|29|l . 



These transition matrix elements and the eigenenergies 
allow one to calculate 6f{E) for arbitrary energies. How- 
ever, the states in Eq. IjBip are offset in energy from 



E 



A(0) 
X 



■ E" . For probe energies E ~ E'^'"' , virtual tran- 



sitions to the states \A+B^), \S~), and iTg") are domi- 
nant, and 9f{E) simplifies to the approximate expression 
given in Eq. (pTT|l . 



\{B+BZ\P-M\'' 
\{A^+'\P+M\' 
\{f,'\P^\eB.+)\' 

\(S-\P-\eB^+)\' 
\{B^-XP-\eB.+)\' 



Pa^b^Ai 
Pb^a^b^ 

2 ^ 

2 ^ 

(l -Ps^A 



(B3a) 
(B3b) 
(B3c) 

(B3d) 
4- (B3e) 



In Table |nl we list all states with two electrons and 
one hole which contribute to Op{E) up to 0{tl). We 
also provide the general expressions for the eigenenergies 
to order 0{t'^) and evaluate them numerically for the 
parameters discussed in Sec. 
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TABLE 11: Eigenstates with two electrons and one hole 
which contribute to the FR angle up to second order in tc- We 
also list the corresponding eigenenergies to 0{t^l) and evaluate 
them for the parameters discussed in Sec. IVl 
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